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ABSTRACT 


A digital computer program for the geometrically nonlin- 
ear static and dynamic response of arbitrarily loaded shells 
of revolution known as SATANS is modified, removing a re- 
quirement that the applied loads be symmetric about a datum 
meridional plane. The ability to include arbitrary 
naturally-occurring initial geometric imperfections of any 
shape, location or extent is added to the program. Two 
specific axially loaded circular cylinders are analyzed 
using actual initial imperfections. Detailed numerical 
studies are conducted to determine the buckling loads of the 
perfect and imperfect cylinders. Several sets of boundary 
conditions are considered. An analysis procedure for imper- 
fect shells is proposed. The results indicate that a means 
"ο conduct a ee analysis of an arbitrarily imperfect 
shell of revolution has been developed. 

A study of the buckling behavior of the clamped shallow 
spherical cap modelled under various boundary conditions at 


the pole is included as an addendum. 
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I. INTRODUCTION 


Shells are used as structural elements in nearly every 
space vehicle, aircraft or other structure where light 
weight is of importance. Consequently, a large number of 
digital computer programs for the structural analysis of 
this shells have been developed. Many of these programs 
were briefly described by Bushnell in Ref. l. An assessment 
concerning the analysis capability based on computational 
systems that have application to a broad class of shell 
structures was presented by Hartung in Ref. 2. In his as- 
sessment, Hartung discusses the limitations, scope, accuracy 
and efficiency of virtually every major shell structural 
analysis program currently available and under development. 
Of particular interest here, is the computer program SATANS 
5. tactic And Transient Analysis, Nom near Shells =- de- 
veloped by Ball for the geometrically nonlinear analysis of 
arbitrarily loaded shells of revolution /Refs. 3 and 4/. 

Pest, in section II, a brief description of SATANS is 
presented including the basic theory, the program capabili- 
ties and limitations, the solution method and procedure and 
some general comments. 

Second, in section III, SATANS-I is constructed. Modi- 
fications to SATANS making the program more user oriented 
are discussed. Included in these modifications are a reor- 
ganization of most of SATANS into a 'storable package, ' 


incorporation of an equipment-independent plotting package 
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which uses a high-speed line printer, modification of itera- 
tion control, discussion of extrapolation technique and 
alteration of boundary load application procedures. 

Third, and most important, the subject of imperfections 
is treated in section IV. Thin isotropic shells of revolu- 
tion have consistently displayed experimental buckling loads 
much lower than the theoretically predicted buckling loads. 
This disagreement has generally been attributed to (1), the 
mathematical inability to correctly model the boundary 
conditions as they exist on the test specimen, and (2), a 
high degree of imperfection sensitivity, i.e., naturally- 
occurring imperfections in the initial shape of the shell 
cause a premature triggering of critical buckling modes and 
early shell collapse. An attempt by Dantone /Ref. 5/ to 
modify SATANS and to include initial imperfections in the 
analysis of cylindrical shells was severely hampered by a 
restriction on the use of SATANS, namely, all applied 
loading and initial conditions -- and hence, initial imper- 
fections -- must be symmetric about some meridional plane. 
Discussion of the removal of this restriction and introduc- 
tion of initial imperfections into the program are covered 
in this section. 

Fourth, and also in section IV, the implementation of 
the initial imperfection capability is presented. The 
resulting program is called SATANS-II. Steps for actual 
Eonstruction of SATANS-II from SATANS-I are presented as 


Ale” instructions for its use. 
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Fifth, in section V, actual imperfect shell analyses 
are performed on two specific axially loaded circular cylin- 
ders. The actual topography of the cylinders has been 
measured by Arbocz and Babcock /Ref. 6/ and this data pro- 
vides the opportunity to include actual naturally-occurring 
initial imperfections into the analysis by SATANS-II. 
Detailed numerical studies are performed to determine the 
classical and consistent bifurcation buckling loads of the 
perfect cylinders and the results are compared with previ- 
ously published results. Stability solutions for the 
imperfect cylinders under six boundary conditions are 
obtained and compared with experimental values. An analysis 
procedure for imperfect shells is proposed based upon these 
results. 

Lastly, as an addendum to this report, the clamped 
shallow spherical cap is investigated under various boundary 
conditions at the pole. While on an industrial tour with 
the Structural Analysis Division of the Palo Alto Research 
Laboratory of the Lockheed Missiles and Space Company, 
Sunnyvale, California, Ryan noted that few shell structural 
analysis programs have the capability of analyzing a shell 
of revolution which has an initial or a final pole -- 
continuous structure through the axis of revolution. This 
investigation attempted to determine the effects on the 
solution caused by substitution of a small hole or a small 


rigid plug for the pole. 
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TES DESCRIPTION OF SATANS 


The computer program known as SATANS is capable of per- 
forming static and transient analyses of thin, linearly- 
elastic shells of revolution subjected to arbitrary loads 
and temperature distributions. It is based on Sanders' 
nonlinear field equations /Ref. 77 for the conditions of 
small strains and moderately small rotations. The program 
can be used under the following conditions: 

1. The geometric and material properties of the shell 
must be axisymmetric, but may vary along the shell meridian. 

2. The shell material must be isotropic, but the 
modulus of elasticity may vary through the thickness. 

3. Poisson's ratio must be constant. 

4. The boundaries of the shell may be closed, fixed, 
elastically restrained or free, and loaded or unloaded. 

5. The applied forces, pressures, temperatures and 
initial conditions must be symmetric about a datum meridional 


plane. 


A.  GEOMETRY AND COORDINATE SYSTEM 

The coordinate system used to locate points within the 
shell is the s, Ө, £ orthogonal system, where s is the 
meridional distance along a reference surface, 0 is the 
circumferential angle measured about the axis of revolution 
and ç is the perpendicular distance from the reference 


surface. The shell geometry and coordinate system are shown 
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in Figure 1. Positive directions for displacements, rota- 


tions, forces, moments and loads are shown in Figure 2. 


Bee METHOD OF SOLUTION 

The governing partial differential equations are reduced 
to an lnfinite number of sets of four second order differen- 
tial equations in the meridional and time coordinates by 
expanding all dependent variables in a sine or cosine series 
in terms of nd, where n is the Fourier index (harmonic) and 
8 is the circumferential coordinate. The sets are uncoupled 
by utilizing appropriate trigonometric identities and by 
treating the nonlinear coupling terms as 'pseudo loads.' 
The meridional derivatives are replaced by conventional 
central finite difference approximations, and the displace- 
ment accelerations are approximated by the implicit Houbolt 
backward differencing scheme /Ref. 37. This leads to sets 
of algebraic equations in terms of the dependent variables 
and the Fourier index or harmonic. At each load or time 
step, an estimate of the solution is obtained by extrapola- 
tion from the solutions at previous load or time steps. The 
sets of algebraic equations are repeatedly solved using 
Potters' form of Gaussian elimination /Ref. 9/. The pseudo 
loads are recomputed after each solution and the procedure 
continues until the solution converges. 

In the static analysis, a load-displacement history is 
obtained by an automatic incrementing of applied loadings 
in equal steps until the number of iterations required to 


achieve convergence in any harmonic exceeds a prescribed 
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limit. The load increment is then reduced by a factor of 
five and the procedure continues until a prescribed number 
of load steps, or a prescribed number of load step reduc- 
tions have been made. 

Since the method of solution is based on the pseudo load 
approach, the shell reacts equally, in a linear fashion, to 
any change in either the applied load or the pseudo load. 
Thus, failure of the solution to converge in any harmonic 
can be attributed to two types of nonlinear behavior, both 
of which are illustrated in Figure 3. Referring to the 
figure, curve A experiences an inflection point or a maximum. 
Either condition causes a zero or nearly zero slope and a 
convergence failure. This curve is typical of softening 
nonlinearity. Curve B, on the other hand, is typical of a 
stiffening nonlinearity and a convergence failure will occur 
when the slope becomes large. It is necessary, then, to 
examine the load-displacement behavior of the shell to 
determine the cause of convergence failure. A detailed 
discussion of the theory and implementation of SATANS can 
be found on pages 8 through 29 and in Appendices B, C and D 
ОТ Ref. 3, and on pages 8 through 28, pages 56 through 61 


and page 100 of Ref. 4. 


foe USER PREPARED SUBROUTINES 

Five small subroutines describing the shell geometry, 
the stiffnesses, the loading (physical and/or thermal), 
and the initial conditions (in the dynamic case) must be 


prepared by the user. Complete descriptions for the 
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preparation of these subroutines can be found on pages 34 


Burough 38 of Ref. 4. 
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LIT CONSTRUCTION OF SATANS=I 


Stein /Ref. 10/ observes that a basic impediment to the 
widespread use of existing computer programs is the difficult 
and laborious process of assimilation of a program's basic 
theory, its implementation and its use by a prospective 
researcher. He also warns the 'unwary analyst' of the dan- 
gers inherent in computer output, which is a 'blizzard of 
numbers.' As a consequence, SATANS is continually being 
modified in an attempt to make it more usable. The most 
recent modifications made for this purpose are described in 


this section. This new version of SATANS is called SATANS-I. 


A. PACKAGING 

SATANS has been reorganized to allow 99% of the code to 
be prepackaged for repeated use by researchers. Compiled 
versions of all subroutines and main controlling routines 
can now be placed on a disk, in a data cell or on a magnetic 
tape, allowing the user to 'submit' only the user written 
subroutines and input data. A complete listing of this 


storable package is given in Appendix E. 


B.  PLOTTING ROUTINE 

A great advantage is gained if the computer output can 
be presented in graphical format. Such plots, however, are 
often generated using mechanical or photo-electronic devices 


which are not universally usable. Consequently, SATANS has 
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been enhanced with an equipment-independent plot package 
utilizing self-contained subroutines and the high-speed 
line printer. 

The main program listed in Ref. 4 could not be modified 
EN ncorporate the logıc and control for a graphical ability 
as it then exceeded the capability of the FORTRAN-H compil- 
er. Therefore, the main program of Ref. 4 has been replaced 
by a new main program and two subroutines described as fol- 
lows: 

:;ΞΙΠ5 1S5 Che new 'main.' It reads all data 
cards, prepares the title page of output and passes control 
to either of two basic controlling subroutines, 

STATIC - This is the controlling subroutine for a static 
analysis. It conducts initialization, establishes geometry, 
stiffness, loading and boundary conditions, controls plot- 
ting of input information and controls the solution proce- 
dure. 

DYNAMC - This is the controlling subroutine for a dynam- 
ic analysis and it operates in a fashion similar to subrou- 
Eine STATIC, 

The graphical ability is obtained by the introduction of 
two controlling subroutines, PLOTl and PLOT2, four plotting 
subroutines, PLOTIT, SCALIT, ROUND and DRAWIT, and modifi- 
cation of subroutine OUTPUT for the selection of output data 


for plotting. A single data card with request information 


1тһе four plotting subroutines are slightly modified 
versions of those contained in Ref. 11 under the listing 
EUSSGUIPLOT." 
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is required to generate the plots. The plotting routines 
will handle a maximum of 100 stations. This routine is 


included in the listing given in Appendix E. 


C. MISCELLANEOUS ALTERATIONS 

The iteration lower limit has been raised to two itera- 
tions. This change was introduced because it was noted 
that immediately after a load step size reduction, several 
consecutive solutions were obtained in one iteration, a 
completely linear process. This was then followed by non- 
convergence and, in some instances, premature program termi- 
nation. Setting a minimum of two iterations automatically 
introduces nonlinearity into each load step and appears to 
Solve this problem, as illustrated in Figure 4. 

A quadratic extrapolation technique was investigated 
using a static analysis of the highly nonlinear shallow 
spherical cap. No gain in solution accuracy was seen and 
a price was paid in program execution time due to approxi- 
mately 10$ greater total iterations required, not to mention 
the increased storage space requirement. This phenomenon is 
exactly as observed by Stricklin et al /Ref. 12/, hence the 
linear extrapolation technique has been retained. 

Since boundary loads are introduced through the boundary 
conditions and since the program considers the boundary 
conditions to be the same in every harmonic, changes indi- 
cated in Appendix B of Ref. 5 are permanently incorporated 


into the program. This alteration causes the applied 
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inc at a boundary to b€ taken in the zeroeth harmonic 
Бу, i.e., when n=0. 

All 'COMMON' statements have been ordered and separated 
and some slightly modified. This alteration allows easier 
reading and correlation of FORTRAN variables with specific 
subroutines. 

The output ability has been expanded to allow printed 
data (and plots) at up to 36 vice six meridians. 

In the case of a shell with an initial (or a final) pole, 
determination of the shear force at the pole in the n=l] har- 
monic has been deemed an inconsequential and lengthy computa- 
tion and has been removed. Appropriate changes have been 


made in subroutines INLPOL, FNLPOL and OUTPUT. 


Pee Wok OF SATANS-I 

Instructions for preparation of input data cards for 
SATANS-I are contained in Appendix C and the initiating main 
program and user written subroutines appear in Appendix D. 
Appendices D and E contain a complete listing of the 
Sompucer program, SATANS-I. 

SATANS-I requires 206 k-bytes of main core when run on 
an IBM-360/67 computer under OS/MVT and compiled by the 
FORTRAN-H compiler at optimization level 2. Execution times 
are the same as for SATANS unless plotting is requested, 
and then is increased approximately 0.01 seconds for each 


100 point plot generated. 
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LVS CONSTRUCTION OF SATANS=I1 


As discussed previously, a serious impediment to the 
consideration of initial imperfections by Dantone /Ref. 57 
was the requirement that the imperfections be symmetric 
about a datum meridian. This restriction results from the 
fact that all dependent variables were expanded in a single 
Sine or cosine series. The analysis presented in this 
section removes this restriction by introducing the use of 
full trigonometric series expansions for the dependent 
variables, and then modifies the resulting new relations 
where required for the introduction of initial imperfections. 
The computer program that incorporates these changes is a 


modification to SATANS-I, called SATANS-II. 


A. INTRODUCTION OF FULL TRIGONOMETRIC SERIES EXPANSIONS 
The original assumption by Ball /Ref. 3/7 was that a 
normal pressure load could be expressed in a Fourier cosine 


series as 


(n) 


q = (0¿h,/E,) pe Соз пө (1) 


Om 8 


where Og the reference stress level, 


Шо (he reference thickness, 

a the reference Young's modulus, 

n = the Fourier index, or harmonic, 
and р = the nondimensional Fourier coefficient of the 
normal pressure load. The normal displacement W 
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corresponding to this pressure was taken as 


W = СО ЫЕ) ү w WU cog nd (2) 
О 


where a is the shell characteristic length and w is the 
nondimensional series coefficient o* the normal displace- 
ment. For a completely arbitrary case, representing q and 


W with the complete Fourier series gives 


p! (B) cos no +) p2 (n) 


1 


q = CA ( Sin nê) (3) 


O wg 


and 


W = (ao /E,) ( ; ы (1) соѕ пө + Y w2(n)sin n0) (4) 
1 


| indicates a previously existing coef- 


where superscript 
wT Et and superscript ^ indicates a new coefficient. 
Similar series expansions can be written for the other 
pressures, displacements, etc. 


The nonlinear terms associated with the full Series are 


memastrated by 67, 
S 


a B (с„/к„)°| ( ) $1 U) cos n6) 24 e $2 0) sin Be 
1 


2 (n) Sin пө) | (5) 


+ 2 ) $1 UJ cos no) ( ) Е 
0 1 ) 


where the Ó < are the nondimensional series coefficients of 
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the meridional rotation, given as 


ο ΕΤ: ἡ 


(6) 


2(n) 


y? (n) = [μ᾽ a 2 


ΙΤ. 
S S 


Originally, in Refs. 3 and 4, the nonlinear term Qe was 


represented as 
С ) me τ ΞΟΠ, ) 81 (n) cos ne (7) 


since 


Cos(n0)Cos(m0) = %Cos(n-m)8 + %Cos(n+m)® 


The 2 term contains the products of all the a's whose 
harmonic sum or difference equals n. There are now two 
additional terms in equation (5). These can be represented 


‚as 


(90/0) ° ( ; $20) sin n6)? ғ о, /Е, 82 (Mcos ne (8) 


and 


ИЕ) фі (2) сов пө) ( ) $2 0) Sin no) 
0 l 


(9) 


IE O) cn n 
1 


12 


where superscript indicates the product of a previously 


existing coefficient and a new coefficient. 
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Thus, 


$2 a ES ( > олсе nd + ) 82 (n) cos no 


(10) 


t) eat Gon no ) 
` 1 


Similarly, the nonlinear term ФЕ (and also ®° since it is of 


the same form) is given by 


92 = Ko T ) De n0) 2+ ( ) $2 cos no)? 


(A) 
+ 2( 0" sin ne) ( ] 02M cos na) | 
1 0 
which can be represented as 
ὃς E ( ) = ng + ) B2 cos no 
(12) 
+ J acum n6) 
1 
The nonlinear term ке 15 given by 
Ф δ, = 52 ( ) фі (п) соз no} { ) Еа n6 
+ (7) ЕЕ no) ( ) ф2 (2) соз пө) 
0 0 
(13) 


+ ( 


$2 sin ne) ( Jed U? sin no) 
1 


-6 18 


$2 U) sin "st 


*í( 


ф2 (0) соз ne) 


-wg 
omg 
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which can be represented as 
© zı(n) y giz(n) 
= о pon 
DM ο δι ) By Sin пө + ) Bug Cos n8 
(14) 
+ 7 82 "sin no) 
O 

Following this scheme, similar expressions can be obtained 
for the remaining nonlinear terms ре PoNor Φ Νορ, ὄρδρς, 
$N, and ΦΝρ. 

All of the nonlinear terms fall into one of the three 
categories Cos(k0)Cos(£40), Sin(k0)Sin(Z80) or Cos(k0)Sin(260). 
Since these three forms were already accounted for in 
SATANS, there was no need to develop any new routines for 
the full series expansions. It was only necessary to intro- 
duce each new nonlinear term into the appropriate existing 


summation in subroutines PHIBET and TEAETA. 


INTRODUCTION OF INITIAL IMPERFECTIONS 
The initial imperfections only influence the nonlinear 
terms in the equations. The strain-displacement relation- 


ships used in this analysis are 


= U” + + (6% + 6° 

ες U W/R. | Š ) 

25 = u A ον s ($° + $°) Gis) 
TE (у^ + U`/r - с/г + ° $.) 
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where 


2 = —W + U/R, 
Pe = -W’/r + V/Rg (16) 
Ф = k(V^ * r^V/r - U'/r) 
and where Ece Egr E 6 represent the reference surface 
s 
strains, Ф, $ $ represent the reference surface rotations, 


and superscript primes and dots represent partial differen- 
tiation with respect to s and 0, respectively. 

Consider a shell whose geometry differs slightly from a 
shell of revolution. For this imperfect shell, a deviation 
from a specific shape occurs in the form of a transverse 
displacement, W. This W produces no strain. Defining W* 


as a total displacement, 


W* = W +W (17) 
and replacing и іп equations (15) with W* of equation (17) 
gives 
ce. + ©. = U « W*/R. * h(os^ + 97) 
S S S ES 
=n ` 2 2 
-- y” ^ JA L (ax 
ο ο ο ενα ив ы ¿(92% + 0”) (18) 
Ze E y + ° = ^ 
cr ;(V^ 4 U'/r - r^V/r + оох) 
where 
"ΠΧ 
oe W + U/R, 
(19) 
k = «ΤΑ 
ФА W* /r + V/R, 
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y Are the pseudo strains due to W only. 


= SU μες ο N, and €. = €4 = Eg = 0. 


and ες, ερ and ες 
ν 5 


When U = 


Substituting these relationships into equation (18), solving 


BOE Ea, ερ, = ss and placing these terms in equations (16) 
lead to 
© * "d ° mi L (W^? 
EU * W*/R. * h(o* + $ TR =) 
(20) 
2 = eo 
Ey = UIT c x UE + W*/R, + д + ^) - Ν/ Χα - (W°) /r 
ερ” "τυ πα ΕΙ ДҮ | 


Expansion of equations (20) making use of equations (16), 


(17) and (19) lead to 


= A i 2 2 B T. 

€g 7 U^ * W/R, * 5(92 + 0%) - {WO} 

EX V'/r-* r^U/r * W/Rg * n(o + Фф?) - (W"0,/r) (21) 
EU ^ ° ^ _ τω, —= 

U + U Zr — r^V/r + 6% - {W°O, + Wo /r}) 


where the terms enclosed in braces are the newly generated 


terms due to the imperfections, and 


ó* = 6 = {W°} 
S S 
(22) 
k = Se 
ο ш (W'/r] 


The presence of imperfections has no effect on the changes 
of curvatures since the curvature-displacement relations are 


linear. 
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Expressing the imperfection partial derivatives -W^ and 
-W'/r as imperfection rotations $, and ο οὓς Strain- 


displacement equations (21) can be expressed as 


= y” ο $ 2 
....υ ην. + {20.9 } + 9°) 


Eg = V'/r * r^U/r + W/R, + EU + 2 ооо 


с (v^ + U /r - Eos * 9 + $ Фф }) 


Ө 50 


where the terms in braces are the newly generated terns. 
Imperfection rotations also affect the nonlinear term, ON. 
(and the others listed after equation (14) ). For these 
Eus the reference surface rotations Ф„ апа Po must be 
replaced by a and A of equations (22). 

To illustrate the procedure which is used in the follow- 
ing section to introduce the initial imperfection rotations, 
consider the single series expansion, equation (7), repeated: 
пеге, 


$2 = ΕΙ à DASS: n0)? z ayy Ea) ce n0 (7) 
0 


Modifying this equation to include the effect of initial 


imperfections indicated by equations (23) results in 


02 + (20,8.) - (oy/E)?| ( ] o1 z 
5 ss HE | ; $.' Cos n6) 


+ 2( Î 4} cos ne) ( J 5! cos ne) (24) 
0 0 J 


_ ο πως 
= S E ) Be © n0 
0 
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A redefinition of В com elude “the cross=product term is 
now required, along with redefinitions of the other f terms. 
and the n terms in the other nonlinear term expressions. A 
redefinition of these terms alone is sufficient for the 
introduction of initial imperfections into a SATANS analysis 
where only the single series expansions are used. This was 


the thrust of the thesis by Dantone. 


C. COMBINATION OF FULL TRIGONOMETRIC SERIES EXPANSIONS AND 
BIITIAFR IMPERFECTIONS 


Consider that for every meridional rotation ὃς, where 


2s = (о/о) } ф1 (2) соѕ ne + ; $2‘) sin no) (25) 


there exists an initial meridional imperfection rotation 


$., defined as 


Ф. = iu BT ) $1 U) cos 'ng + } 92 N sin пө) (26) 
0 — 


The nonlinear term Ф, previously given by equation (5), now 
becomes (62 * 26.6.) , according to equation (23), where o* 

is given by equation (5) and 945,9, is the product of equations 
(25) and (26), as illustrated for the single series expan- 
Sion in equation (24). All the nonlinear terms of equations 
(15) have been converted to those in equation (23) for the 
full series expansions and the resulting relations are 

given in Appendix A. “he definitions of the B terms (given 


previously in Ref. 3) have been modified to include the 
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effects of both full series expansions and initial imperfec- 
tion presence. 

The remaining nonlinear terms (those listed after equa- 
tion (14) ) also require modification but are of simpler 
form since only the first power of the rotation is involved. 
A redefinition of the n terms to include the effects of full 
series expansions and appropriate initial imperfection rota- 
tions has also been made. 

In this manner, the validity of the strain-displacement 
equations (equations (5) of Ref. 3 and equations (15) of 
this report) and the equilibrium equations (equations (8) 
of Ref. 3) is retained. Introduction of full series 
expansions and initial imperfections in the nonlinear terms 
was accomplished by a redefinition of the 8 and n terms, 
which are the nondimensional coefficients for the nonlinear 
terms in the strain-displacement and equilibrium equations, 


respectively. 


D. SUBROUTINE MODIFICATIONS 

Introduction of the full series into the program was 
implemented by use of negative index numbers for the new 
sine or cosine portions of the expansions (the terms with 
the superscript ?). With the requirement that the axisym- 
metric harmonic (n=0) be included, the program was expanded 
to handle a maximum of 25 modes, i.e., the zeroeth harmonic 
(n=0) plus twelve other harmonics (positive and negative). 
Each positive index and each negative index will be referred 


to as a mode, i.e., the harmonic n=5 has two modes, -5 
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and +5. A typical correspondence between harmonic, mode 


and Fourier index is given in Table I. Harmonic zero (n=0) 


and any combination of up to twelve harmonics may be 


included in the 
Table I are the 
negative values 
the full series 
the first mode, 


first mode of a 


analysis. The positive values of NN in 
modes analyzed in SATANS and SATANS-I. The 
of NN are those modes added to incorporate 
expansions in SATANS-II. Note that n-0 is 
and that the negative value is always the 


pair of modes. This is a requirement. 


TABLE I 


HARMONIC, MODE AND FOURIER INDEX CORRESPONDENCE 


HARMONIC| MODE| FOURIER| FORTRAN ASSIGNMENT 
а INDEX STATEMENT IN SUBROUTINE 
PLOAD OR TLOAD 
0 1 0 


NN (1)=0 





-1 NN(2)2- 
NN(3)=+1 
NN(4)=-2 
NN (5) =+2 
NN (6) =-5 
NN (7) =+5 
NN (24) =-18 
NN (25) =+18 
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Seven subroutines of SATANS-I required changes to 
incorporate the full trigonometric series expansions and 
the initial imperfections. 

Subroutine PHIBET calculates the $'s and carries out 
the summation procedures for calculating the P nonlinear 
terms. Since extensive modification was required, the 
subroutine was altered such that those sections handling 
the full trigonometric expansions and imperfections are 
entered only if required, i.e., if a full series analysis 
1s necessary. 

Subroutine TEAETA calculates the inplane forces and 
the n nonlinear terms. It was also extensively modified, 
again, utilizing new sections only if required. 

The subroutine MODES, a highly complex procedure for 
determining modal coupling required minimal alteration 
since each unique mode number is stored at every, or every 
other storage location, depending on use of SATANS-I or 
SATANS-II, respectively. Each unique harmonic is accessed, 
then, by an appropriate step size change in the controlling 
loops. 

The subroutine OUTPUT required the addition of a section 
to include the new half of the trigonometric series in the 
generation of a summed solution at a specified meridian. 

The subroutines POLE, INLPOL and FNLPOL compute the 
solution at an initial or a final pole. Due to the compu- 
tation of pole conditions based on n, which now has positive 


and negative values, these routines were modified 
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accordingly. Effectively, the only change in the equations 
presented in Appendix D of Ref. 3 involve a replacement of 
n by -n where appropriate. 

A complete listing of these seven modified subroutines 


is given in Appendix F. 


fee SATANS-II 

To modify SATANS-I to handle completely arbitrary 
loading and initial imperfections, the following steps 
must be performed. 

l. Replace subroutines PHIBET, TEAETA, MODES, OUTPUT, 
POLE, INLPOL and FNLPOL with those in Appendix F. 

2. Replace the affected COMMON statements in Appendices 
D and E with those listed in Appendix F. 

3. Add subroutine IMPERF, also listed in Appendix F. 
This subroutine contains the values of the imperfection 
rotation coefficients, $. and Pa- Instructions for actual 
preparation of these coefficients will be discussed in de- 
tail in section V. If imperfections are not to be consider- 
ed, the Ф. апа Фе ате never laceessed. 

4. In Appendix E, make the following additions after 
the card indicated: 


Subroutine After Card Number Insert 


DYNAMC 00002270 CALL IMPERF 

STATIC 00002230 CALL IMPERF 

PMATRX 00001100 TEN UN). 19, 0)P(2,1,19)=-1%0 
PMATRX 00002000 IF(N(MN).LT.0)2FP1(2,2)-1.0 
SATANS 00000380 READ(5,108)JUMP,MPERFS 
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SATANS 00001400 108 FORMAT (215) 
Sen Appendix E, modify the following cards: 


Subroutine Card Number Should Read 


DYNAMC 00000800 DIMENSION VBAR(25) ,AVB(25) 
PMATRX 00000900 NN=IABS (N (MN) ) 
PMATRX 00001640 NN=IABS (N (MN) ) 


6. In Appendix C, modify the data preparation instruc- 
tions by adding a data card after card #4 according to the 
following format: 


Card Column Format Item Interprétation 


4.1 1-5 Т5 JUMP 1 = single series expansion 
2 = double me expansion 
4.1 6-10 I5 MPERFS 0 = no imperfections 
1 = imperfections 


An initial imperfection solution requires the use of full 
series in the analysis. Specification of the parameter 
JUMP indicates whether the solution will be as in SATANS-I 
or will utilize the full trigonometric series expansions of 
SATANS-II. Specification of the parameter MPERFS indicates 
whether the analysis will include initial imperfections. 
ДЕЕР спав if MPERFS equals 1, JUMP must equal 2, but if 
JUMP equals 2, MPERFS may be either 1 or 0. 

7. Add imperfection data just prior to the last card of 
the data deck. This data is to be placed in arrays PHIXB 


End PHITB. 
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F. USE OF SATANS-II 

After completing the changes in section IV.E, SATANS-II 
exists as a program capable of performing all the functions 
of SATANS-I with the following restrictions: 

1. The number of meridional stations cannot exceed 
100, (KMAX < 101). 

2. The number of modes cannot exceed 25, (MNMAX < 26). 

3. The product of meridional stations and modes cannot 
exceed 951, (KMAX*MNMAX « 952). 

4. Plotting ability is not available for the imperfec- 
tion input data. 

Specification of harmonics to be considered (accomplish- 
ed in subroutine PLOAD or TLOAD) must include the axisym- 
metric harmonic (n=0) and both negative and positive entries 
for any other harmonic, as indicated in Table I. 

The instructions for the preparation of the input data 
are given in Appendix C.  SATANS-II requires 426 k-bytes of 
main core on an IBM-360/67 computer when run under OS/MVT 
and when compiled by the FORTRAN-H compiler at optimization 
level 2. Execution times vary widely from four minutes for 
a 35 station and 3 mode analysis, up to 30 - 60 minutes for 


a 35 station and 25 mode analysis. 


G. VALIDITY OF SATANS-II 

SATANS-II was checked for its ability to obtain a cor- 
rect solution using the full series expansions. This check 
was made using a clamped cylindrical shell subjected to 


pressure loadings as described in Figure 5. A solution was 
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C. ASYMMETRIC LOAD 


&- = SOLUTION COMPARISON POINTS 
SATANS-IL CHECK LOADINGS 
FIGURE 5 
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obtained for the loading shown in Figure 5A with SATANS-I 
and the summed solution was printed at meridian Ө = 0°. The 
cylinder was then analyzed by SATANS-II for the following 
three orientations: 

1. Loading defined as in Figure 5A - a symmetric load. 

2. Loading defined as in Figure 5B - an anti-symmetric 
load. In this case, the summed solution was printed at 
meridian 8 = 90°. 

3. Loading defined as in Figure 5C - a general load. 
In this case, the summed solution was printed at meridian 
Е 459. 

Load definitions for the four check runs appear in 


Table II. In each case, the solution obtained by SATANS-II 


TABLE II 


FOURIER COEFFICIENTS FOR CHECK LOADINGS 


FOURIER COEFFICIENTS OF APPLIED LOADING 





ΜΟΌΕ |------------------------------------------- 
SYMMETRIC ANTI-SYMMETRIC NONSYMMETRIC 
0 B 

-1 0.0 -45.015 

+1 -63.66 : 452015 

-3 0.0 EIS 065 

+3 '- : -15.005 

-5 ϱ.0 ΓΡ» 

+5 -12.73 : 792003 

-7 0550 -6.431 

+7 292095 š +6.431 
Note: the coefficients are based on p. 7 50 psi. 
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agreed with the SATANS-I solution ehroughir rour significant 


digits. 


“This reference to solution includes the load-deflection 
history, the ultimate nonconvergence level, and the displace- 
ment, force and moment distributions. 
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V. IMPERFECT SHELL ANALYSIS 


А. PREPARATION OF THE IMPERFECTION DATA 

With the developments of section IV implemented, 
SATANS-II requires only the Fourier coefficients of the 
imperfection data to conduct a nonlinear analysis. Recall 
Eon equation (26) that these Fourier coefficients, $. and 
$6 are determined by knowledge of the actual imperfection 
Fotations, Du mange К МКМ кес атас from equation (22) that 


9 


the imperfection rotations are given as 


Q = -{7 
s 


(29) 
$. - -W'/r 


0 

This means that the imperfection rotation Fourier coeffi- 
cients can be determined by knowledge of the initial normal 
Mertlection deviation, W. This imperfection information can 
be obtained from a mathematical description, a graphical 
plot or from actual experimental measurement. A simplified 
data generation flow chart is presented in Figure 6. 

1. Data Generation by Mathematical Description 

In a note on imperfect shallow spherical cap buck- 

ling, Kao /Ref. 13/7 analyzes a cap subjected to a mathemat- 
ically defined dimple imperfection existing over one-quarter 


of the shell. The dimple is given by 
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GRAPHICAL EXPERIMENTAL 


POT MEASUREMENT 





MATHEMATICAL 
COMPUTATION 


=| 





i Kn Taf) n) eU 


SIMPLIFIED IMPERFECTION DATA GENERATION 
SCHELE 


FIGURE 6 
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W = 64(s/L}*{i=s/L)*hécos 26,  -459«0« 459 


(28) 
we 0, A60025 


where L is the total length of the shell meridian and where 
s/L ranges from 0 at the pole to l at the edge, and ô is an 
intensity factor, 0 <= ó «Ξ 1. The geometry of this problem 
is shown in Figure 7. Kao used a two-dimensional finite 
difference mesh of eight meridional and 16 circumferential 
stations in his analysis. A numerical model of, for example, 
35 meridional stations and 25 modes would yield a SATANS-II 
solution of considerable accuracy. The discrete initial 
N deflection deviations W at each mesh point need not 
be computed, as partial differentiations of equation (28) 
give expressions directly solvable for discrete values of 

Ds and DE OE CEN 6S and nper ection rotations at loca- 
tions exterior to the affected area (-459 < 9 < 459) аге 
zero. To illustrate, the meridional imperfection rotations 
that exist along the line ABCDEA of Figure 7 would appear 

as curve A of Figure 8. The discrete values generating 

this curve can easily be expanded into a full Fourier series 


(n) inserted into SATANS-II 


and the resulting coefficients $_ 
for the analysis. 

The imperfection of Figure 7 is symmetric about the 
meridian 0 = 0°. As an example of a more complex imperfec- 
tion, another dimple of the same form could be located in 


the region defined by 1009 « 8 < 1909. In this asymmetric 


case, the meridional imperfection rotations along line 
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8=45° 





PN 
0-190? 1 
m. 
0= 325° 
SHALLOW SPHERICAL CAP 


IMPERFECT QUARTER WITH 
DIMPLE IMPERFECTION GIVEN BY 


W/hz 64x" (1-XP^ d Cos 20 
WHERE X=st 


MATHEMATICAL IMPERFECTION GENERATION 
FIGURE 7? 
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A.SIMPLE DIMPLE MERIDIONAL IMPERFECTIC 4 ROTATION 


B С D E A 


1/2 3/2 


B. COMPLEX DIMPLE MERIDIONAL IMPERFECTION ROTATION 


MERIDIONAL IMPERFECTION ROTATION 
REPRESENTATION 


FIGURE 8 
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ABCDEA of Figure 7 would now appear as curve B of Figure 8. 
The discrete values generating this curve can also be easily 
expanded into the full trigonometric series and its coeffi- 
cients inserted into SATANS-II. 
2. Data Generation by Graphical Plots 

Graphical representations of W can be easily con- 
structed and the patterns are near limitless in design or 
location. Theoretical or actual creases, furrows or dents 
of any shape need only be graphed, the discrete imperfection 
displacements determined, the discrete imperfection rota- 
tions DE and er computed via a finite differencing tech- 
nique, the Fourier expansions performed and the resulting 
coefficients inserted into SATANS-II. 

3. Data Generation by Experimental Measurement 

Arbocz and Babcock /Ref. 6/ have developed a means 
of mapping the topography of circular cylindrical shells, 
and can thus describe the deviations of an actual shell. 
The procedure, described in detail in Ref. 6, produces 
imperfection data in the form of nondimensional normal 
displacements from a perfectly circular cylinder. For the 
sample shells considered here, 49 data points were taken 
around each circumference and 29 equally-spaced circumfer- 
ential traces were made. The first and last trace were 
0.64 inches from the ends of the shell and the spacing 
between traces was 0.24 inches.  Pertinent geometrical and 
material data concerning the shells of interest here is 


ESHmnEarned in Table III. 
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TABLE III 


IMPERFECT SHELL OVERALL CHARACTERISTICS 


SHELL DESIGNATION] лез А-14 


Radius 4.00" . 00" 
‘Length ΝΟ. 
Young's Modulus x 10*, E |15.10 psi |15.80 psi _ 
Wall thickness, h | 0.00444" | 0.00437" — 
Number of axial traces | 29 | 29 . 
‘Axial increment, A Î 0.24" | oze . 
“Number of data points | 1421 | | 1420 —— 





TO prepare the imperfection data on shell A-14 for 
use by SATANS-II, the nondimensional deviations were first 
dimensionalized by the wall thickness. To retain the 
original axial increment of 0.24 inches between circumfer- 
ential traces and to make the data compatible with a 35- 
station meridian, the shell length was increased by 0.08 
inches at each end, or a total of 2% (See Figure 9). This 
required three additional traces at each end of the shell 
corresponding to stations l, 2 and 3, and 33, 34 and 35. 
Data for the additional six circumferential traces were 
generated via the modified linear extrapolation procedure 
illustrated in Figure 10. At each of the 1,715 data points, 


the meridional and circumferential rotations were computed 


52 














v 
Ф £ 
£ O 
a E 
О (ο 
O Е 
Е 27 3 e 
28 | 31 
SHELL A-14 
5 à 23 
34 
35 > 
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MERIDIONAL STATION NUMBER CORRESPONDENCE 


FIGURE 9 
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e - KNOWN DATA POINT (REF.6) 

E TWO STATION AVERAGE 

A =EXTRAPOLATED DATA POINT 

Wij= IMPERFECTION DISPLACEMENT 
OF THE ith STATION AND THE 
jth MERIDIAN 


NOTE: NUMBERS CORRESPOND TO 
SATANS-IL MERIDIONAL STATIONS 


IMPERFECTION - FREE MERIDIAN 
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DATA POINT EXTRAPOLATION TECHNIQUE FOR W 


FIGURE IO 
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using central finite differencing in the appropriate direc- 
tion. After nondimensionalization according to the 
SATANS-II scheme, each circumferential vector of 48 
meridional and 48 circumferential imperfection rotations 

was expanded via a periodogram into a Fourier series accord- 


ing to 


—(0) 25 —c (n) 285 
(Е0/00)Ф. = p £ ) Ó < Cos n0 + ) 


—s (n) 
1 ὃς 


Sin no (20) 


and 


(E¿/0,) 9, = E 4 3 $$ cos nO + b $5) Sin nos) 
E his point, the Fourier coefficients 6. and $6 are of 
proper form for direct insertion into SATANS-II. The 
computer program which performed the actual imperfection 
data processing for shell A-14 is given in Appendix G. 
4. Selection of Important Harmonics 
Since the number of meridional stations was 35, the 


program maximum of 25 modes can be used in the analysis.” 


Note, however, that in equations (29) and (30), there are 
49 modes describing the imperfections. Therefore, a selec- 
tion of a maximum of 25 imperfection modes must be made. 
The harmonics to be considered in the analysis should in- 
Mame those Critical to the shell and those significant 


in the imperfection. The description of the harmonic 


3Recall the restrictions given in section IV.F. 
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Selection process for shell A-14 is given in a later 


Section. 


B. BOUNDARY CONDITIONS 

The testing configuration used by Arbocz and Babcock 
1 axially loaded cylinder is shown in Figure ll. The 
ends of the shell are attached to an end ring and a load 
cell with cerrelow, a low melting-point alloy. The ends of 


the shell are believed to be clamped so that 


W = $ =V = 0 (31) 


The axial load on the test specimen was applied by means of 
four screws located at the corners of the square end plates. 
At each load increment, the screws were adjusted to give as 
orm a circumferential load distribution as possible. 

In all versions of SATANS, the boundary conditions are 


represented by the matrix equation 


N U 
N. V 

ay, 208007 = {4} (32) 
Q W 
η 
ὃς | Msi 


Boundary conditions can be specified on N. or U, N "ΓΟ 


50 
Qs, or ana Ф. ог М. by prescribing input values for the 
matrices / Q / and / A /. There is one equation (32) for 


Sach boundary of the shell. In their present form, SATANS-I 
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and SATANS-II assume the column matrix (2) as zero for all 
harmonics except the axisymmetric harmonic. 

As was mentioned earlier, boundary condition discrepancy 
between mathematical model and physical experiment is cited 
as a major contributor to thin shell buckling load disagree- 
ment. For this reason, three interpretations of Figure 11 
were made. The cylinder was represented as experiencing an 
axisymmetrically applied force Νε on poth ends (condition A), 
an axisymmetrically applied force N. on one end and a 
meridional deflection restraint on U at the other (condi- 
tion B), and an axisymmetric incremental end shortening U 
at one end with meridional deflection U restrained at the 


4 These three conditions were investi- 


EUST (condition C). 
gated on the clamped cylinder, and on a simply supported 
cylinder for further comparisons and validity checks. 


Mathematical assignments for boundary conditions A, B and 


C are shown in Table IV. 


C. ANALYSIS 
T. Buckling of the Thin Perfectly-Circular Cylinder 
A necessary part of this attempt to develop a pro- 
wor predicting the buckling load of imperfect shells is 
the establishment of the program's ability to predict the 
buckling load of perfect shells. Consequently, this section 


compares the results of a SATANS-II analysis of the thin 


4 ; . "E 
It 1s believed that boundary condition C best repre- 
sents the actual boundary conditions in Figure 11. 
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cylinder without imperfections with the bifurcation buckling 
analysis presented by Almroth /Ref. 14/. Almroth's study 
deals with the influence of edge constraints on the critical 
load of cylindrical shells of the type used by Arbocz and 
Babcock. Equation (33) is given by Almroth for the axisym- 
metric displacement shape of the simply supported cylinder 


Meeundary condition A) prior to buckling. 


W = r(2vyN. E EMI ansia tas) Sanitas) 
(33) 
+ A,Cos(a,s)Cosh(a,s)) | 


where 


а = «(1 S N.) š 


-h 


κ -(2γ) 


h/r/(12(1 - v*))* 


. 
t 


A 
«4 
| 


s =-N./[2yEnh) 


το, (El) 


rO | 
I 
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(11 - E2) 5Sin(o,)Sinh(a;) - N Cos (a,) Cosh(a3)) 


A, = - E — o oc 
(1 - N2) '5(Cosn? (a,) + Sin? (a,)) 
B ( (1 - N2}®cos(a,)Cosh(a,) 2 N Sin(a,)Sinh(a;)) 
= -= 
— n р  — — Ə—əƏ, s,zsqp 
[el E N*) * (Cosh* (o2) - Sin? (a,)) 
Ме (a4,LD)/ (2r) 
BN (a2L)/ (2r) 


and where No is the axial load, p is a uniform normal pres- 
Smee, fh is the thickness, E is Young's Modulus, v is 
Son's ratio, L ls the length and r is the radius. This 
solution is referred to as the consistent prebuckling solu- 
tion, as opposed to the membrane or classical buckling solu- 
tion where W is a constant along the meridian. Equation 
(33) has been solved for shell A-14 and W is plotted as a 
MENCEdOn Of s in Figure 12 for an Es = 0.21904. Also shown 
in Figure 12 is a plot of W versus s obtained from SATANS-II 
for the same МИ for both 35 and 69 meridional station anal- 
yses. Note in Figure 12 that even though only a few sta- 
tions lie in the boundary zone, the trend is towards the 
rect solution. The displacement of the 35 station solu- 
tion, however, is only a rough approximation of the actual 
Skate. 


TGS ND OAS determined Ehe minimum bifurcation 


buckling load, Мест, and the corresponding value for n for a 
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range of parameters L/r and r/h, using both the consistent 
prebuckling displacement W given by equation (33) and the 
membrane or classical prebuckling solution, where W is 
constant, in the Donnell equations. For the simply support- 
ed cylinder with an r/h of the same order of magnitude as 
πι A-14 and an L/r = 0.5, the critical load is N. „0.863 


(N =39.5 lb/in) with n=27 for the consistent analysis, and 


ser 


Neo roi. 0 (чы 746 lb/in) with n=27 for the membrane pre- 


buckling analysis. Thus, the consideration of the axisym- 
metric edge effect, which is analogous to an imperfection, 
reduces the buckling load by 14.1%. Parametric studies 
conducted by Almroth show usage of from 100 to 800 merid- 
ional stations over the half-length of the shell. 

As shown in Ref. 5, use of the membrane prebuckling 
displacement in Sanders' equations for the simply supported 
cylinder simplifies the problem to an eigenvalue problem. 


The buckling mode is 


W = 6 Sin(mts/L) Cos (n8) (34) 


The equations for the matrix elements in the eigenvalue 
problem are given in Appendix B. A plot of the critical 
values of m and n is given in Figure 13. The minimum buck- 
ling load occurs when n=27 and m=15, and corresponds to an 


N mm) 25>) 1 b/in. 


Scr | 
SATANS ability to predict the minimum bifurcation 
load of a shell for a specified value of n has been demon- 


Strated by Stilwell /Ref. 15/7. SATANS-II analyses of shell 
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A-14 (simple support, boundary condition A) were conducted 
over a wide range of n for both the consistent and membrane 
prebuckle cases. Figures 14 and 15 give plots of the crit- 
ical buckling load versus the number of meridional stations 
for various values of n. The dashed parts of the curves 

are estimated extrapolations. It is clear that as the 
number of stations is increased, SATANS' critical buckling 
load decreases and the associated harmonic increases. While 
the n=27 curve predicts the highest buckling load over the 
range of stations used in both cases, extrapolation of 

these curves indicates that the n=27 curve would predict 

the lowest buckling load if the number of stations exceeded 
200 to 300. Note that the membrane prebuckle critical 
buckling curves of Figure 15 seem to flatten much more 
quickly than for the consistent case. This appears to 
confirm the occurence of a reduction in critical buckling 
load due to the presence of an axisymmetric edge restraint, 
Eei ed by Almroth. Figures 16 and 17, which present plots 
of the critical buckling load versus critical buckling mode 
of shell A-14 for the 35, 69 and 100 station analyses, also 
predict this lowering of buckling load due to edge restraint. 
the 'fanning' of curves in Figure 16 indicate an increase 

in the number of meridional stations will produce a lesser 
buckling load decrease than in Figure 17. These figures 
also predict a lower critical buckling load and a higher 


critical buckling mode with increasing meridional stations. 
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Figures 14 through 17 indicate that, given enough 

M ions, SATANS-II can predict a critical buckling load 
and mode for the thin-walled shell that agrees with 
Almroth's prediction. However, for the 35-station analysis, 
which is the analysis to be used for the imperfect shells, 
ehe critical load of 51.1 1b/in and harmonic of 16 are 
Significantly different from the correct values. Asa 
result of the coarse mesh used by SATANS-II in the following 
analysis, the critical buckling loads for the imperfect 
shells will be higher than the experimental values. 
Nevertheless, it will be shown that the introduction of a 
totally arbitrary, naturally-occurring imperfection pattern 
on the imperfection-free mathematical model of shell A-14 
= zules in a significant reduction in critical buckling 
load. This is not only sufficient to demonstrate the 
capability of the program to perform the analyses, but it 
also provides an estimate of the degree of imperfection 
sensitivity of the shell. It will be shown that application 
of this estimate of sensitivity to Almroth's buckling load 
provides an accurate prediction of the buckling load of the 
imperfect shell. 
2. Dre, of Imperfections 

when prepared as indicated in section V.A.2, the 
imperfection rotation Fourier coefficients are in a form for 
direct introduction into SATANS-II. Figure 18 shows the 
applicability of the periodogram method of expansion of the 


imperfection data according to equations (29) and (30). 
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Curve A is the graph of a solution of equation (30) pre- 
senting the circumferential imperfection rotation at 
station 10 plotted against 0. It is seen to pass smoothly 
through all the data points. 

Only the zeroeth harmonic (n=0) and twelve others 
can be accommodated by SATANS-II. In the case where the 
shell response and the imperfection are limited to a few 
major harmonics, all important harmonics can be included in 
the analysis. For example, the ۸=6 shallow spherical cap 
has only the n=2 harmonic as the major harmonic of response, 
and the imperfection curves of Figure 8 require only a few 
Fourier harmonics for adequately accurate representation. 
Also, the few required harmonics for imperfection represen- 
tation are in the vicinity of the n=2 harmonic. This pro- 
blem is will within the limits of SATANS-II. 

The thin-walled cylinder, however, responds in a 
large number of harmonics. This can be seen in Figure 16 
-- many harmonics produce critical buckling loads close to 
the value of the minimum critical buckling load.  Addition- 
ally, the imperfection requires many harmonics. Thus, as 
mentioned in section V.A.4, a difficult selection of harmon- 
ics for use in the analysis must be made. Six harmonics 
were selected on the basis of the shell response and six 
harmonics were selected on the basis of the imperfection 
characteristics. Of course, all 13 selected harmonics will 
contain both the imperfection characteristics and the shell 


response characteristics. 


m3 





For the shell response harmonic selection, the 
critical buckling loads of shell A-14 for specified values 
of n, discussed in the previous section, are presented in 
Table V. Values were assigned to each harmonic according 
to the legend of Table V and six were selected based on 
these assignments. 

For the imperfection harmonic selection, the sums of 
the absolute values of the Fourier coefficients for both 
sine and cosine terms of both meridional and circumferential 
imperfection rotations were tabulated and are presented in 
Table VI. Values were again assigned to each harmonic ac- 
cording to the legend of Table V, and six additional harmon- 
ics were selected. 

The complete harmonic selection process for shell 
A-14 was accomplished by a comparison of the values assigned 
to each harmonic for both the imperfection and the shell. 
Table VII presents this comparison along with the actual 
harmonics selected for the analysis. Note that no overlap 
of important harmonics occurred, as the imperfection re- 
quired harmonics seven and below and the shell required 
harmonics ten and above. The harmonics used with shell A-13 
were n-0 and n-2,3,4,6,8,9 for the imperfection and n-10,11, 
12,13,14,15 for the shell response with a simple support, 
and n-13,14,15,16,17,18 for the shell response with a fixed 
support. 

Once the specific harmonics for inclusion in the 


analysis are determined, the appropriate ὁ, and %, 
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TABLE V 


CRITICAL BUCKLING LOADS OF INDIVIDUAL HARMONICS 


БИЙКЕ О SURPORT OTMELL SUBPORT 


HARMONIC BUCKLING|P/Par ASSIGNED|BUCKLING Dar | ASSIGNED 
LOAD | VALUE LOAD VALUE 
x Pg 


E —— 


neq ee 


11| 


w CE | κο — 2 





_ 23_ 1.340 _ 1.290 
24 1.344 
Notes: 1l pert 50 р51. 
2. The trigger in each harmonic is an internal 


pucccupo cu US psi: 


3. Each analysis included the indicated harmonic 


plus n=0. 
4. Value assignments are according to the following 
legend: SYMBOL | VALUE 


absolutely necessary 
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TABLE VI 


RELATIVE IMPERFECTION ROTATION COEFFICIENT SIZES, SHELL A-14 


Φρ 
B pun; 
VALUE| SIN |VALUE| COS |VALUE| SIN |VALUE 


о OO YA UA UNO 


Notes: 





х 
* 
+ 
+ 
+ 
* 
+ 
+ 
+ 


Numbers represent the sums of the absolute values 
of the imperfection rotation coefficients over 
the entire shell. 


Value assignments were made according to the 
legend of Table V. 
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TABLE VII 


HARMONIC SELECTION, SHELL A-14 


SHELL SELECTED 
RESPONSE HARMONICS 









ETXED 
SUPPORT 





SIMPLE 
SUPPORT 


FlXED SEMELE 
SUFEORT | SUPPORT 


PRIMARY PRIMARY 
ШКОЛ Мн ee NATE | 
X X X 
X X X 
к | ж X X X 
* | + X X X 
+ | * X X X 
+ | * x X X 
T ΞΡ 
* | + X X X 
+ = 
£o - X 
= = * = X Р 
+ | + * + X X 
- | + κά * X X X 
- | - X X X 
- | - * * X X X 
- * κ κ x X 
- = * * X X 
= = * * X 
= zs * 
= T 
= + 


Harmonic zero is required by SATANS-II. 
Value assignments are taken from Tables V and VI. 


ΙΧ’ indicates a selected harmonic. 
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coefficients (obtained from the periodogram) must be loaded 
into the SATANS-II arrays PHIXB and PHITB (performed in 
subroutine IMPERF) in the following order: 

l. The first KMAX locations contain the $.'s (or фо '5) 
of mode 0. 

2. The second KMAX locations contain the AE (or $o's) 
of mode -i. 

B. The third KMAX locations contain the $.'s (or $o 's) 
of mode +i. 

4. The fourth KMAX locations contain the $.'s (or φο' 5) 
of mode -j. 

P οετετα πει all modes are entered. 
KMAX is the FORTRAN name for the number of meridional sta- 
tions to be included in the analysis. One special case must 
be considered where the shell has an initial or a final 
pole. In this case, rotations are zero in all but the п=1 
harmonic. If this harmonic is included in the analysis, 
the imperfection rotation coefficients at the first (or 
last) station should all be set to zero except for the 
coefficients in the n=l harmonic. In the case of an initial 
pole, the coefficient associated with the -l index should be 
placed in the first location of the second set of KMAX loca- 
tions in arrays PHIXB and PHITB, and the coefficient assoc- 
lated with the +1 index should be placed in the first 
location of the third set of KMAX locations. For a final 
pole, the coefficient associated with the -l index should be 


ШО е in the first location of the next to last set of KMAX 
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locations, and the coefficient associated with the +1 index 
ΠΠ be placed in the first location of the last set of 
KMAX locations. 

To illustrate the effect of the selection of har- 
' monics on the imperfections, equation (30) was solved using 
only the 13 selected harmonics. The resulting circumferen- 
tial imperfection rotation representation at station 10 
appears as curve B in Figure 18. To determine the effect of 
the inability to include all nine harmonics deemed important 
to the shell response with a fixed support, an alternate 
selection of harmonics was made for this case (see Table 
VII). A complete SATANS-II analysis was performed with 
this alternate selection and the results were found to 
differ from those using the primary harmonic selection by 
approximately 2%. 

3. Results 

' Six analyses (boundary conditions A, B and C with 
clamped and simple supports) were made for each shell where 
the imperfections were reduced by a factor of 1000. These 
'minute' imperfections acted only as triggers for the modal 
deflections and the resulting solutions gave the shell's 
consistent critical bifurcation buckling loads. These 
results are presented in Table VIII under the category of 
minute imperfections. Selected cases were analyzed where 
E minute imperfection triggers were replaced by a 10 'psi 


* 


internal pressure trigger and the results were identical. 
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TABLE VIII 


DUCK ENG TOADS, MINUTE VERSUS, FULL, IMPERFECTIONS 











SHELL A-14 SHELL А-13 


REDUC- |BUCKLING | REDUC- 
TION LOAD TION 
(%) ΓΙ» Τη) (©) 


EWBPORT| IMPERFEC- | BOUNDARY |BUCKLING 
CONDITION| LOAD 


(1Ь/1п) 


С 

А 

М 

Р 

Р — 

D Full 19055 
ευ 

5 Minute 

Т 

М 

Р 

L 26.30 

E Full 24.90 





* insufficient data. 


The same six analyses for each shell were also made 
using the full imperfections. Results for the final load at 
which a converged solution could be obtained are also shown 
in Table VIII along with the percent reduction in the buckling 
load for each case. 

Since the full imperfection is introduced at the 
first load step as a pseudo load, the program requires con- 


siderably more iterations than normal to obtain the first 
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solution. The subsequent solution usually proceeds in a 
completely normal fashion. However, convergence difficul- 
ties may be encountered in cases where the imperfection 
level is high. This was the case with shell A-13, as some 
of its imperfection deviations exceeded 3% shell wall thick- 
nesses, and fluctuated over 6 thicknesses in a 8 arc of only 
539. For shell A-13, a 70% initial imperfection level was 
the maximum for which a converged solution could be obtained 
on the first load step. Consequently, buckling loads were 
computed with a maximum of 80%, 85%, 90% and 95% of the full 
imperfection data. For the 80% imperfection, 70% was intro- 
duced at load step one and 10% was added at load step two. 
For the three other solutions, 5$ increments were added at 
various load steps. Since SATANS uses the total value of 
the variables and not the incremental values, each inter- 
mediate solution is correct for its given load and imper- 
fection level, but the load-deflection characteristics are 
invalid. An imperfection level of 95$ was the absolute 
maximum for which a converged solution could be obtained for 
for shell A-13, and a short 5$ extrapolation was made to ob- 
tain the data presented in Table VIII for this πο. This 
procedure is illustrated for the shell with a simple support 
in Figure 19. A comparison of load steps and iterations for 
both shells under these conditions is presented in Table IX. 
For shell A-13 with clamped boundary condition A, a con- 
verged solution was not obtained for an imperfection level 
above 85%. This condition is therefore not included in the 


results. A similar problem was encountered by Dantone Гвен 
5/ with this boundary condition. 


Sal: 
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TABLE IX 


LOAD STEPS VERSUS ITERATIONS 


SHELL A-14 SHELL A-13 






LOAD|INCREMENTAL |ITERATIONS|LOAD|INCREMENTAL |ITERATIONS 
SEF IMPERFECTION | REQUIRED STEP IEEE CLON TV REOUTRED 






























ADDITION ADDITION 
1 100% 7 1 70% 24 
2 2 2 10% 6 
3 2 3 5 16 
4 3 4 12 
5 3 5 9 
6 3 6 3 
7 3 7 8 
8 3 8 6 
9 3 9 5% 19 
10 3 10 3 
L1 3 11 5% 60 
12 4 12 58 
13 4 13 18 
14 4 14 16 
15 5 15 16 
16 7 16 12 
L7 10 
18 16 1S load step size reduction 
15% load step size reduction| 17 | 4 
19 3 2nd load step size reduction 
20 2 
18 3 
2nd load step size reduction| 19 2 
21 | 3 34 load step size reduction 
3rd load step size reduction termination 
22 | 2 (Dee 10.6 δη τη) 
termination 
CF 46.0 Ib/in) 
Note: This data corresponds to boundary condition B, simple 


EUPDPOREL, 
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bee DISCUSSION 

The SATANS-II investigation of the buckling load of the 
perfect cylinder (discussed in section V.C.1) indicates that 
a numerical analysis using 35 stations results in a higher 
puedrcted buckling load than the correct analytical buckling 
load. Consequently, the SATANS-II results for the predic- 
tion of the critical buckling load of the imperfect cylinder 
are expected to be high. Table X compares the results of 
the SATANS-II analysis (column A) with the actual experi- 
mental results of Ref. 6 (column E) and shows this to be 
the case. The use of more stations to lower the predicted 
buckling load can be accomplished only by either decreasing 
the number of harmonics considered in the analysis or by 
increasing overall program size. Both of these choices are 
impractical at the present time since the number of harmon- 
ics already in the analysis appears to be low considering 
the number required for shell response and the imperfection 
data, and the present program core storage requirement is 
large. 

It is possible, however, to apply the SATANS-II results 
for the percent reduction in the buckling load due to the 
imperfections (see Table VIII and Column B, Table X) to a 
more accurate buckling load of the perfect shell from an 
analysis like that of Almroth to get a more accurate estimate 
of the buckling load of the imperfect shell. When the per- 
cent reduction from the consistent minute imperfection 


buckling load to the consistent full imperfection buckling 
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load, determined by SATANS-II, ls applied to the accurate 
Consistent bifurcation buckling load of Almroth (column C, 
Table X), the result is a prediction of the actual buckling 
load of the imperfect shell. These imperfect shell predict- 
ed buckling loads are shown in column D of Table X and are 
seen to be very close to the actual buckling loads of the 
imperfect shells, given in Column E of Table X. 

The percent reductions in buckling loads determined 
by SATANS-II are considered to be applicable to the accurate 
consistent bifurcation buckling loads of Almroth since the 
imperfection characteristics are seen to be concentrated in 
the lower ten harmonics (see Table VI), where shell response 
does not appear to be largely affected by the number of 


meridional stations in the analysis (see Figures 14 and 15). 


E PROPOSED ANALYSIS PROCEDURE 

When the SATANS-11 predicted buckling load of the perfect 
shell is higher than the analytical prediction, as is the 
case with the extremely thin-walled cylinder, the following 
procedure is proposed. 

1. Determine the consistent bifurcation buckling load 
as close as possible. 

2. Determine the buckling load with SATANS-II using 
minute imperfections. 

3. Determine the buckling load with SATANS-II using 
full imperfections and compute the percent reduction from 
step 2. 


4. Apply this reduction to the buckling load of step 1. 
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This procedure is expected to work only when the result 


of step 2 is reasonably close to that of step 1. 
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VI. CONCLUSIONS 


A. SATANS-I 

The computer program SATANS-I has been constructed from 
SATANS. It is a more user-oriented version of the program 
by Ball, whose primary modification is the introduction of 


an equipment-independent graphical output ability. 


B. SATANS-LI 

SATANS-II is a further modification of SATANS where the 
I riction that all applied loading and initial conditions 
be symmetrical about some meridional plane has been removed. 
In addition, SATANS-II possesses the capability of including 
experimentally-determined naturally-occurring imperfections 
or malformations induced during fabrication, mathematically 
described imperfections of any type, frequency or location, 
or graphically described imperfections due to damage suffer- 
ed after fabrication. The imperfections are introduced into 
SATANS-II by using the imperfection deviations along discrete 
circumferences to determine the associated meridional and 
circumferential imperfection rotations at the mesh points, 
performing a Fourier expansion on these rotations and then 
use of the Fourier coefficients. In cases where the impor- 
tant shell response harmonics and the harmonics required for 
imperfection representation exceed the limits of SATANS-II, 


selection must be made. 
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Imstructions for the use of both SATANS-I and SATANS-II 
are contained in Appendix C and complete listings of both 
programs are given in other appendices. SATANS-I requires 
210 k-bytes (53k words) of core on an IBM-360/67 computer 
using 'single precision,' while SATANS-II requires 426 
k-bytes (107k words) of core. Imperfect shell analyses 
performed by SATANS-II in this report used a 35 station 
meridian and 25 circumferential modes. Execution times 
were an average of 60 minutes, although execution times 
of between 15 and 30 minutes can be expected for analyses 


using 15 or less modes. 


C; ШАРЕКЕЕСТ SHELLS 

SATANS-II has demonstrated an ability to perform a 
numerical structural analysis on the arbitrarily imperfect 
shell of revolution. Imperfections of over three shell wall 
thicknesses in depth have been handled. Increasing solution 
difficulty relates directly to increasing imperfection 
complexity and increasing shell instability. In those cases 
where extremely thin-walled shells result in perfect (imper- 
fection-free) buckling loads from SATANS that are higher 
than analytically determined buckling loads, a procedure has 
been devised where the percent reduction due to the intro- 
duction of full imperfections can be applied to the analyti- 
Cally determined buckling load, and an imperfect shell 
buckling load predicted. These predictions have been shown 
to be very close to the actual imperfect shell buckling 


loads for the two specific axially loaded cylinders considered. 
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ADDENDUM 
PEDI GITAL COMPUTER STUDY OF THE EFFECTS OF VARIOUS ARTIFICIAL 
BOUNDARY CONDITIONS AT THE POLE OF A SHALLOW SPHERICAL SHELL 


1% INTRODUCTION 


The design of many shell of revolution structures re- 
quires continuity of the structure at one or both poles. 

At such a pole, the radius becomes zero and special 
conditions must be imposed on the deflections and the stress 
resultants to assure finite strains and equilibrium. Asa 
consequence of the difficulty of incorporating these 
conditions into a structural analysis computer program, 
many researchers model such a structure as having a small 
hole or a small rigid plug at the pole. This addendum 
investigates the applicability of the various artificial 
pole conditions by examining their effect on the stability 
of a clamped spherical cap subjected to axisymmetric and 
asymmetric pressure loads. 

Various polar artificialities are discussed and the 
boundary equations associated with their use are given, as 
are those for the continuous pole. A determination of the 
applicability of the various poles is made by comparisons 
of their effect on the shell buckling loads and by comparing 
the errors in the normal displacement and meridional mem- 


brane force distributions incurred by their use. 
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The computer program SATANS, developed by Ball /Ref. 3/7, 
was used to conduct this study and selected results of the 
SATANS analyses are checked for validity with the computer 
programs BOSOR /Ref. 15/ and STAGS /Ref. 15/. SATANS has a 
special routine that incorporates the proper conditions to 
be enforced at a pole.  BOSOR also has this special ability 
but STAGS is not yet so ο... 

SATANS, a one-dimensional finite difference equilibrium 
formulation based on Sanders' equations, was run on an 
IBM-360/67 computer under OS/MVT.  BOSOR, a one-dimensional 
finite difference energy formulation based on the Donnell 
equations, was run on a UNIVAC-1108 computer under EXEC-VIII, 
while STAGS, a two-dimensional finite difference energy 


formulation was run on a CDC-6600 computer under COMPASS-12. 


Cas er bebruary, 1972, 
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IT: STATEMENT OEF THE PROBLEM 


A. GEOMETRY AND LOADING 
The structure investigated is the clamped shallow 
spherical cap whose geometry can be specified by the single 


nondimensional parameter à, where 


A = 2(3(1-v2))* (n/n) * (AD-1) 


аша where H is the rise of the shell and h is its thickness. 


For this study, Poisson's ratio was taken as 0.3, Young's 


ΙΙ. as 30 x 10° psi and h as 1.0 inches. Shells of 
À = 6 and A = 18 were investigated however only the results 
of the À = 6 shell are given in detail. 


The classical buckling pressure of a complete sphere is 


given as dg where 
a ^ 
ο wale ο) (AD-2) 


and where Βς 15 the radius of curvature. Fox this study, 
two loading conditions were used. The cap was subjected to 
a uniform normal pressure (symmetrical load) and a uniform 
normal pressure over half the cap (asymmetrical load). 
These loadings are shown in Figure 20. In each case, the 
magnitude of the applied pressure is gq and the buckling 


results are presented as the ratio σα. 
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B. POLE CONDITIONS 

A shell is considered to have a pole if the structure is 
continuous through its axis of revolution. At a pole, the 
radius is zero and special conditions must be imposed on 
deflections u, v and w, and on the meridional moment m. to 
assure finite strains and equilibrium. The computer program 
SATANS is equipped with the capability to impose these 
special boundary conditions. A full derivation of the 
required conditions is given in Appendix D of Ref. 3. The 
results of that derivation are presented below. 


For the zeroeth Fourier harmonic (n=0): 


u = = = = 
V Mon t0 0 
MN m. ^ = 6) 
(AD-3) 
m = m 
sd © 
tg 7 tg — b(1-v) (u^tuswe58..) -t., 
For the first Fourier harmonic (nz1l): 
(utv) = w = u” = 0 
Mo = m. = Mag = 0 (AD-4) 
ts Е Eg ~ Ese =" 


When a double sign (+) appears, the upper sign applies for 


an initial pole and the lower sign applies for a final pole. 
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От the second harmonic (n=2): 


ay = w= mM = 0 
L et 
Ө Ө 
° å (AD-5) 
My = ms = ims 
ο δι 
For all higher harmonics (n>2): 
u = у = W = 0 
πο τ m. = 50 = 0 (AD-6) 
ο Ec = Ese = 0 


Bor rotations at the pole, for the first harmonic (n=1): 
Ф = -w tW uU 
$6 = +P (AD-7) 
> ov su j 
and for all other harmonics, 
Cameo DO 
š Е (AD-8) 
ф = 1%(2v” + nu ) 


Numerous artificial boundary conditions can be substi- 


tuted for the correct conditions at the pole. For example, 
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a small hole can be placed at the pole and can exist with 
the boundary completely free, or have conditions imposed 
such that the meridional rotation bs and deflection u are 
restrained. A meridional line load ο ος magnitude suffi- 
cient to replace that 'lost' due to the hole can be applied. 
Another artificial pole condition is that best described as 
a small rigid plug. Here, both the meridional as well as 
the circumferential deflection is restrained along with the 
meridional rotation. A normal line load -P whose total 
equals that loading 'lost' due to the hole or plug can be 
applied at the boundary of the hole or plug. 

Preliminary numerical studies showed that the line load- 
ings effected a noticeable change in the solution only in 
the immediate area of the boundary, hence, the line loads 
were not considered further. Preliminary studies also 
showed that the cap behavior exhibited only minor changes 
when the small hole with a partially restrained edge was 
substituted for that with a totally free edge. This 
investigation, therefore, concentrated on two artificial 
boundary conditions at a pole and on the continuous pole. 
The artificial pole conditions used and their associated 
mathematical expressions are: 


l. The totally free hole: 


Ns = 0 ; Noe = 0 ; Qs = 0 ; Ms = 0 (AD-9) 


U ο ο ο το, >. = 0 (AD-10) 





C. TYPES OF ANALYSIS 

Three types of buckling analyses were conducted on the 
cap under each of the three pole conditions and these were: 

1. The nonlinear axisymmetric analysis of the symmet- 
rically loaded clamped spherical cap. 

2. The bifurcation buckling analysis of the symmetri- 
cally loaded clamped spherical cap with nonlinear prebuckling 
axisymmetric displacements. 

3. The nonlinear collapse analysis of the asymmetrically 
loaded clamped spherical cap. 

In the artificial pole cases, hole and plug radii were 
taken as 0.17%, 1.70% and 3.30% of the shell meridian. The 
SATANS analyses used a 100 station meridian for the axisym- 
metric cases, a 66 station meridian with three harmonics for 
the bifurcation buckling cases and a 40 station meridian 


with five harmonics for the nonsymmetric collapse cases. 
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TITS DISCUSSION 


The results of the SATANS analyses are presented in 
Table XI for all cases. Also presented are the normal 
deflections of a point located at 47.5% of the meridian 


which is near the maximum displacement for each case. 


A. THE X = 6 SHALLOW CAP 
ie rhe Axisymmetric Analysis 

In Table XI, the buckling load is seen ES be basi- 
cally unaffected by the altered pole conditions. The normal 
displacement, however, is seen to show a marked sensitivity 
to the pole conditions. Figure 21 shows this sensitivity. 
In the figure, the upper curve is the deflection shape of 
the cap at incipient buckling with a η... pole, and 
the lower curves represent the percent error induced in 
the normal displacement at incipient buckling by the intro- 
Eun of artificial poles. The 0.17% plug and the 0.173 
hole behaved identically. Holes, in general, produced an 
Mage 20% error in the normal displacement and this error 
existed over the entire shell. The plugs, preserving the 
continuity of the structure, produced an average 10% error, 
but again, this error extended over the entire shell. 
Figure 22 shows the behavior of the meridional force Ns for 
' Οοπειπιοας5 pole case, and the N, error due to the 
ШОО icial pole conditions, also at incipient buckling. 


Large errors existed in all cases over the first third of 
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TABLE XI 


- 


SHALLOW CAP (A=6) BUCKLING LOADS AND DISPLACEMENTS AT A POINT 





ANALYSIS MERIDIONAL NORMAL 
TYPE CONDITION | STATIONS / q/q,  |DEFLECTION 
HARMONICS AT 47.5% 
MERIDIAN 
_ و‎ {inches 
| | - 0.939502 _ 
Pole 307: 1.0610 Be 5720 
| 100/1 8. 0200 1.3200 
Axisymmetric |---------- pam - - - - - - -------- w 
Buckling -0.17% 100/1 0.9720 1.1600 
- Hole-1.70% 100/1 0 5690 1.0560 
Symmetric -3.30% 100/1 0.9620 1.0290 
O A A A a == 
-0.17 100/1 0.9720 172590 
Plug-1.70$ 100/1 0.9770 1.1630 
-3.305 10071 0.9900 1.2290 
= O27 750* = 
Pole 3073 0:7851 0.6386 
Go 3 (00752 0.6188 
Bifurcation  |----------4----------4--------t---------- 
Buckling -0.17$ 66/3 0 7119 0.6356 
_ Hole-1.70% 66/3 0.7835 0.6562 
Symmetric -3.30% 66/3 0.7392 0.6264 
PAC 7 aaa 
-0.17% 66/3 0.7719 0.6356 
Plug-1.70% 66/3 0.7726 0.6438 
-3.30% 66/3 027713 0.6513 
_ O57 2007; = 
Pole 40/5 0.5391 0.9830 
Te ي ي‎ O 
Analysis ο μυς 40/5 035363 1.0780 
- Hole-1.70% 40/5 0.5389 1.0060 
Asymmetric -3.30% 40/5 0.5962 1.4110 
Шз Шы сс — C 1 — س س ي س ص‎ O SS SS ص‎ 
-0.17% 40/5 0.5978 0.9478 
Plug-1.70% 40/5 0.6290 0.8100 
-3.30% 40/5 0.7576 0.7248 


m Huang /Ref. 16/ 


t Famili and Archer /Ref. 17/ 
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the meridian and were an average of 10% or less for the 
remainder of the meridian. 
2. Bifurcation Buckling Analysis 

In this analysis, the shell buckling load appears to 
be insensitive to artificial pole conditions for the case of 
the plugs, and shows some sensitivity in the case of the 
holes. This is surprising since bifurcation buckling of the 
À = 6 cap occurs in the Cos(20) mode. While equations (AD-8) 
show that τ, ΟΙ ΕΠΙ ο Ele pole in this mode, rotas 
tions certainly do exist at 1.7% and 3.3% meridian due to 
the buckling mode. The equations for the rigid plug's 
boundary conditions (equations (AD-10) ) specifically pre- 
vent this rotation, while those for the hole (equations 
(AD-9) ) do not. 

Examination of Figures 23 and 24 show that the dis- 
placement and force distributions of the plugs at incipient 
buckling to be less erroneous, again, than the poles. In 
' Figure 23, an average 10% displacement error occurs over the 
entire shell for the plugs, while holes of 1.7% and 3.3% 
meridian cause large errors over much of the shell. This 
trend is repeated in Figure 24, although the errors are not 
as large. From this analysis, the indications are that the 
1.7% hole is a poor artificial pole and a 3.3$ hole is 
unacceptable as it causes a 4.6% error in buckling load and 
over 40$ error in some regions in both normal displacement 
and meridional force distributions. Larger hole sizes would 


apparently increase these errors. 
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ERE NoOusymmetric Collapse Analysis 


The results for the collapse analysis are also pre- 
sented in Table XI, and in Figures 25 and 26 for the 0 = 0° 
meridian, which is symmetrically located under the loaded 
Eon. In Table XI, it is shown that the buckling loads of 
the cap with rigid plugs are quite erroneous (10.85, 16.65 
апа 41.4% for the 0.17$, 1.7$ and 3.3$ plugs, respectively) 
while those for the small holes are quite reasonable. In 
Figures 25 and 26, it is seen that the small holes cause an 
average error of 10% in both the normal displacement and 
meridional force distributions at incipient buckling over 
the entire shell, while the errors caused by the rigid 
Bugs is large. 

In an asymmetrically loaded shell, large rotations 
will exist at the pole and the restrictions imposed by the 
rigid plug appear to make it an unacceptable artificial pole 
for use in this case. 

Figures 25 and 26 also indicate that a 1.7$ hole is 
the best artificial pole for this case. Assuming such a 
'best' hole size exists for another given shell or loading 
condition, an investigation similar to that presented here 
would have to be performed to find it. This is an inherent- 
ly expensive and time consuming process and there is no 
guarantee that there is a 'best' hole size without its 
radius becoming so small that mathematical difficulties 


arise during the solution process. Consequently, the use of 
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artificial pole conditions in nonsymmetric problems must be 


approached with caution. 


B. THE A = 18 SHALLOW CAP 

A deeper shell with A = 18 was investigated in a similar 
fashion and produced similar results. The shell buckling 
load was negligibly affected (except in those cases already 
mentioned for the A = 6 shell) but displacement errors of up 
to 30% and meridional membrane force errors of an average 
10-15% were noted. In general, the à = 18 shell appeared 
Sensitive to the type of artificial pole (except in the 
BEuxudetric case), but as sensitive to any artificial pole as 


me shallower à = 6 shell. 


COMPARISONS WITH OTHER RESEARCHERS 

me Indicated in the introduction, various SATANS results 
were checked for validity with BOSOR and STAGS. These com- 
parisons are presented with others in Table XII and very 


favorable agreement is seen to exist. 
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Iy: CONCLUSIONS 


eron buckling cases if the radius Or the artificiality is 
kept at or below 33% of the meridian. In the case of the 
asymmetrically loaded Shell, this insensitivity only applies 
when small holes are substituted for пе Bole. η general, 
the following additional Conclusions can pe drawn. 

l. Sizable errors are introduced into the normal dis- 
Placement and meridional membrane force distributions by the 
introduction of an artificial pole. 

КО The displacement and force error distributions are 
sensitive to the size and type of artificial pole and are 
usually greater than 10% with some Cases of over 30%, and 


they exist Over the entire shell meridian. 


fOr Symmetrically loaded shells, Apparently their 
continuity-preserving Characteristics greatly outweigh their 
E C Onal-restricting Characteristics, 

4. Small holes appear to be the best polar models for 
asymmetrically loaded ΕΙ. Here, the rotational restric- 
tions of the rigid Plug make its use unacceptable. πο Έρις 
type of loading, an investigation will be required to deter- 


mine the Optimum hole Size, 
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5. In general, the greater the displacement, or the 
higher the membrane force that would exist at a pole, the 
greater the error introduced by the use of an artificial 


pole. 
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APPENDIX A 


NONLINEAR TERMS 


The nonlinear terms of equations (15) are given below. 
Terms involving the initial imperfection rotation terms 


along have been omitted as they cause no straining. 
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1 0 


* 2( > $! UJ sin no) ( Y Ф2 (9) сов ne) | (A-3) 
0 


ni 


E) в1 (9) іп пө + y g? cog пө) 
1 0 


11 3 





фі (1) соз 


a) 

ср 
г 
— 


фі r. 


D 

ср 
== 
mo 


92 MW sin n8) ( 


ر2 

ср 
ر‎ 
ὑπ 


ΝΕ; 


rJ 
CD 
EN 


Ico 
1م‎ moe пө) ( 
ф2 2) 5іп п) (у 5, sin ne) (А-4) 


ф2 60) 5р n9) { 


Сос 


(m) 
Е 


е, 

CD 
--- 
omes 


$1 UJ cos ne)( 


ΝΤ ηθ) { 


$2 0 sin no) ( 


ел 7-78 өмл 8 et-38 -r-38 ‚гл o ο 18 - 8 eoig oig ec-38 
ot-18 -t-318 0618 2-58 20678 r618 et-38g -wmm 8 omg -oig omg -mmg 


= 0 /E, ( ) 81 cos пө +) ἔτ nê) 
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APPENDIX В 
EIGENVALUE PROBLEM MATRIX ELEMENT DEFINITIONS 


The elements of the /A / and /B / matrices used in the 


eigenvalue solution are defined below. 


E Gi, - ETO 
A12 F29: 


mos — F,,0 


Азу 7 -F5,,0 

СЕ а 6. (Be) 

А, = G5, t EOS T G u ο 

Es πμ; 

A32 = G33 7 DESC 2 Gn T Foo An 

Ass = Gss 7 το Т (G,, š Be ο. ü ο ο πο, 
where 

Q = mT/L 


and where 


Bi, = n*/4 

B,. = no/4 

Em = 0 

B,, = n*/4 

Ba, = 02/4 (8-2) 
B,, = 0 

ДЕ 20 

в,, = 0 

B,, = 0? 


CT 





and where all elements of the 4 x 4 matrices / E 7, / F 7 


and / G / are zero except for the following, defined as 


Е = b 
(l JD (Еа (855) 
πίε 
521 (В=5) 
E33 = 2(1:ν) ἀπ' /Βρ 
E34, = 1 
Ει = -d 
Bae 3(l-V)bn/R, = 3(1-v)dn/(8R,) 
F3 = bv/R, - 5 (1-v) n*/Rà 
(B-4) 
Dog = -P12 
F31 = -F33 
Cii = (Eon Ες - (1-v)dn^/(8Rj) 
EX bn /R, - (lv) (Itv)An’/R, 
G23 = -bn/Ro - (1-v) (1+v)an*/Ra 
s = νη/Ερ 
Bee — G23 (B=5) 
G33 = Eh 
Gyo = vdn/Ró 
Gua = vdn*/R$ 
Guy = -1 
and where 
b = {(1-v) (1+v)}7! 
(B= 6) 
ον = (b/12) (h/r) * 
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APPENDIX C 


INPUT DATA GUIDE FOR SATANS-I AND SATANS-II 


Card Columns Format Item Example Meaning 


1 


i 


5 


16-20 


r 25 


Ὁ 


1 5 


36-40 


41-45 


18А4 


Bp 


I5 


I5 


I5 


E5 


15 


I5 


БЕ — Enter any 72 characters 
NO I The problem number, 
0<NO<10000. 
SDYNAMC F For a static analysis, 
set SDYNAMC = F. 
T For a dynamic analysis, 
set SDYNAMC - T. 
IMODE 0 For no modal output data. 
L For modal output data. 
NDIMEN 0 Dimensional output data. 
1 Nondimensional output. 
NTHMAX 8 Summed solution will be 
printed at NTHMAX merid- 
ians, O«zNTHMAX«-36. 
IFREO 3 Solution will be printed 
at the first station, 
every subsequent IFREO 
station and the last sta- 
tion, O«IFREQ«-KMAX. 
ITPRINT 2 Every IPRINT converged 
solution will be print- 
ed. 
IBCINL -l If the shell has a pole 
at the first station. 
0 If the shell has no pole 
aC Chie «lest station. 
IBCFNL -1 If the shell has a pole 
als the last station: 
0 If the shell has no pole 


al ле last station. 


ТЕУ 





Gard Columns Format Item Example Meaning 


2 46-50 I5 KMAX 35 Number of meridional sta- 
tlons. 
Note: KMAX<201 for SATANS 
-I without plots and 
КМАХ<101 for SATANS-I 
τετ or for SATANS 
-II. 


2 51-55 15 MNMAX 7 Number of series coeffi- 
clents used to describe 
the initial conditions, 
pressure and thermal 
loads (and initial imper- 
fections if using SATANS 
SEE): 

Notes FOr SATANS=I: 
O<MAXM<11 and 
КМАХ*МАХМ<201. 

For SATANS-II: 
O<MAXM<26 and 


KMAX*MAXM<951. 
2 61-65 I5 LSMAX 1 For a linear analysis. 
ЭЭ Use many load steps for 
a nonlinear static analy- 
sis. 
750 For a dynamic analysis, 


LSMAX is the number of 
time increments, where 
LSMAX-T,. ,/AT. 


2 66-70 Js LCHMAX 3 The number of load step 
size reductions in a 
static analysis, recom- 
mended range = 2-4. 


0 For a dymanic analysis. 
2 71-75 I5 ITRMAX 1 For a linear analysis. 
30 The number of iterations 


at a load or time step. 
For a nonlinear analysis, 
suggested range = 10-30, 
up to 50 for special 
cases. 


2 76-80 15 EC 0 Initial conditions. ¿Set 
COROS а Static añaly= 
S15, OI for a dynamic 
analysis where the shell 
is at rest at T=0, 
l For a dynamic analysis 
with initial conditions. 
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mara Columns 


3 


3 


ΙΙ; 


12-24 


24-36 


3-48 


49-60 


61-72 


Format Item 


ET2.3 


pu S 


2.5 


F lo) 


οι 


NU 


SIGO 


ELAST 


TKN 


CHAR 


TEEO 


Example 


0.3 


1000.0 


1,0 


2318 


1.0 


Meaning 
POISSON'S ratio, v. 


Reference stress level, 
σο. 

τα νο input data is 
dimensional. 


Reference modulus of 
l sir TI Το. 

ΠΕ ἔπε input data 15 
dimensional. 


Reference thickness, h 
It the input data 15 
dimensional. 


Characteristic shell 
dimension, a. 
Ti the input data is 
dimensional. 
таба static analysis. 


Reference time, Το. 


gun a Q Qm Q a] m m dimus dinum w W q w Q r sm W  --  IIƏ:?) h w A es N N w w q Q q s Q N w w q q A s Q a s m m — «πο πο «πω 


4 


13-24 


212.3 


DELOAD 


EPS 


1235-6 


ο ο 


For a static analysis, 
DELOAD is the load incre- 
ment. It remains un- 
changed until thé solu- 
tion fails to converge 
in ITRMAX iterations, 
when it is reduced by a 
factor of five. A maxi- 
mum of LCHMAX such re- 
Juel ions αλ οσους. 

For a dyhamic analysis, 
DELOAD is the nondimen- 
sional time increment. 


The convergence criterion 
recommended range of 
ОКО < 065 0.001. 


—— — аниц Q Q af a w 7 w q w Q DU E s مس دس‎ a w w s W ee A Å O e A A A A e A جس‎ απ ου «ας دس‎ m. Ñ s s s 


Card 4A is only included for a SATANS-II analysis. 


4A 


d 5 


15 


JUMP 


1 


2 


14.9 


For an analysis using 
single series expansions. 
For an analysis using 
double series expansions. 





era Column Format Item Example Meaning 


4A 5-10 I5 MPERFS 0 An analysis without im- 
perfections. 
1 An analysis with imperfec- 


tions. 
Note: if JUMP=2, MPERFS 
may be 0 or l. 

if JUMP-1, MPERFS 
must be O0 

if MPERFS-1, JUMP 
must be 2. 


Include as many cards 5 as necessary to specify NTHMAX 
meridians. If NTHMAX equals 0, omit card 5. 

5 172 0512.3 "ο. ο A list of circumferential 
coordinates 0, in degrees 
and tenths, where the 
solution printout is de- 
sired. The list must have 


NTHMAX entries. 


If IBCINL--1, omit cards 6 through 14. ТЕ IBCFNL=-1, omit 
wS 15 through 23. Cards 6 through 23 describe the 
boundary conditions at the first, and then at the last sta- 
tion. The boundary conditions exist on the total variables, 
not on the individual harmonics.  Loadings applied through 
specification of boundary conditions are taken in the zeroeth 
EB rmonic (n=0) only, as the column matrix {2%} is set to zero 
for harmonics greater than zero. The boundary conditions are 
EXuensoional. The format of cards 6 through 23 is 4E16.8. 





ι 
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05,15» Card 7,16 Card 8,17 Card 9, 18 
241,1) О ЛЕ?) ΠΡ, 5) Q(1,4) Ng 
22,1) 22) ME Q(2,4) Мв E 
mos, 1) 1005025) (2091.3) Ω(3,4) OQ 
Q(4,1) 00422) 0 (4,3) 0 (4,4) ФЕ 





EE οἱ Card 11,20 Сага 12,21 Сага 13,22 


A(1,3) 
16253) 
163,53) 


ιά 


Card 14,23 
A(1,4) U & (1) 
A (2,4) V 2 (2) 
A(3,4) W аз) 
л (4,4) M. @ (М) 


1, 1) 1 02) 
21) πι; 2) 
5, 1) A(3,2) 
О πι 
Card 24 is: 


Included for a SATANS-I static analysis. 
Included but blank for a SATANS-I dynamic analysis. 
Omitted for a SATANS-II analysis. 


fara Column Format Item 


24 


24 


For the remainder of card 24 entries, 


j| 2 


L2 


L2 


SPLOTS 


SMODAL 


Example 


F 


T 


Meaning 


Indicates plots are NOT 
desired. 

Indicates plots are 
desired. 


Indicates plots are for 
summed solutions only 
Indicates plots are for 
modal solutions only. 


0 indicates that no 


plots are desired for the particular item, and 1 indicates 


that they are desired. 
аасаєеа item versus the station number. 
insure IFREQ-l1l. 


is desired, 


24 


24 


24 


24 


24 


24 


Ὁ Ὁ 


2 


I2 


12 


I2 


T2 


I2 


IRADII 


] GAMMA 


IOMEGS 


IOMEGT 


IDEOMS 


ΤΕ 
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All graphs are plotted as the in- 


If a complete plot 


Plot the radii as computed 
by subroutine GEOM. 


Plot p°/p as computed by 
subroutine GEOM. 


Plot 0, as computed by 
subroutine GEOM. 


Piloto- as Computed by 
τος τς СЕОМ. 


Plot w, as computed by 
subroutine GEOM. 


Plot the stıffness d as 
computed by subroutine 
BDB. 





Card Column Format Item 


24 


24 


24 


24 


24 


24 
24 


24 


24 


24 


24 


24 


24 


24 


24 


17-10 


19-20 


21=22 


23-24 


25-26 


27-28 


Zo 30 
I 32 


33-34 


95-56 


67-358 


39-40 


41-42 


43-44 


45-46 


47-48 


I2 


I2 


T2 


T2 


T? 


I2 


I2 


το 


I2 


I 


I2 


12 


I2 


T? 


12 


I2 


ΙΡ τα 


IBBSTF 


IDDSTF 


TER 


TPS 


ТРТ 


TT 
IMT 


ΕΕ 


IDMT 


INS 


TITI 


INSTH 


IOS 


IMS 


IMTH 


Example 


1 


22 


Meaning 


Plot the stiffness d as 
computed by the subroutine 
BDB. 


Plot the stiffness db/dZ 
as computed by subroutine 
BDB. 


Plot the stiffness dd/dg 
as computed by subroutine 
BDB 9 


Plot the normal component 
of the pressure load. 


Plot the meridional com- 
ponent of the pressure 
load. 


Plot the circumferential 
Component of the pressure 
lead, 


Plot the thermal load. 
Plot the thermal moment. 


Plot d/dz of the thermal 
load. 


Plot d/dz of the thermal 
moment. 


Plot the meridional mem- 
brane force distribution. 


Plot the circumferential 
membrane force distribu- 
Eion 


Plot the meridio- 
circumferential membrane 
πο ο S τορυιςιος. 


Plot the transverse force 
distribution. 


Plot the meridional moment 
distribution. 


Plot the circumferential 
moment distribution. 





аата Column Format 


24 


24 


24 


24 


24 


24 


24 


49-50 


51 5 


5354 


55-56 


DOS 


ου 


61-62 


I2 


I2 


12 


12 


12 


T2 


12 


Item 


IMSTH 


IU 


IV 


IW 


TPHIS 


IPHIT 


TENT 


Example 


Meaning 

Plot the meridio- 
circumferential moment 
dISEEIDULION_ 


Plot the meridional dis- 
placement distribution. 


Plot the circumferential 
displacement distribution. 


Plot the normal displace- 
ment distribution. 


Plot the meridional rota- 
pIonvduastribution. 


Plot the circumferential 


Ἱ ατιο distribution. 


Plot the meridio- 
circumferential rotation 
απο τομ. 


Insert imperfection data here for a SATANS-II analysis with 


imperfections. 


Insure format of the imperfection data is 


compatible with that specified in the user-written sub- 
routine IMPERF. 


23 


Indicates this is the 
only run. 

Indicates another run is 
to be made. Add another 
complete set of data 
cards after this card is 
IRNAGN-1. 
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